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Abstract 

It is shown that the VT symmetric Hamiltonian with the periodic potential V{x) = Acos'^x + 
4iVosin2x can be mapped into a Hermitian Hamiltonian for Vq < .5, by a similarity transformation. 
It is also shown that there exist a second critical point of the potential V{x), apart from the known 
critical point Vq = 0.5, for Vq ~ .888437 after which no part of the eigenvalues and the band 
structure remains real. Relevant physical consequence of this finding has been pointed out. 
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Quite recently, the prospect of realizing complex VT symmetric potentials within the 
framework of optics has been suggested |ll-llOl|. What makes this possible is the formal 



equivalence between the quantum mechanical Schrodinger equation and the paraxial equa- 
tion of diffraction. In this optical analogy the complex refractive index distribution plays 
the role of the optical potential. The parity-time condition implies that the real index profile 
should be even in the transverse direction while the loss or gain distribution must be odd. 
Thus complex optical VT potentials can be realized by judiciously integrating an asym- 
metric gain or loss profile on an even index distribution In these studies it was shown 
that optical VT symmetric materials can lead to behavior that is impossible in standard 
systems. Such effects include double refraction, power oscillation, eigenfunction unfolding, 
non-reciprocal diffraction patterns, bloch oscillations 0] to mention a few. Recently 

Guo et al [Sj] have demonstrated experimentally passive VT symmetry breaking within the 
realm of optics. In ref. the first observation of the behavior of a VT optical coupled 



system that involves a complex index potential is reported. Energy band formation 
and spontaneous symmetry breaking thresholds have been investigated in a few VT sym- 
metric models even in presence of nonlinearity [l| and lattice disorder jsl. In particular, 
while discussing energy bands due to a complex periodic potential, the authors of ref.(ll| 
argued from the graphs of A(i?), i? G M corresponding to these potentials, that some energy 
bands would appear and disappear under perturbations. Later it was shown that the 
appearance and disappearance of such energy bands imply existence of non real spectra due 
to the complex deformation of real intervals. 

In this short note we shall concentrate on the complex periodic potential 

V{x) = Acos^x + 4iVosin2x (0.1) 

where V{x + 71") = V{x). This potential is invariant under the dual action of the space 



reflection operator V and the time reversal operator T [15|. In general the action of V is 
defined by the relation p — —p, x — )■ —x ( p, x stand for momentum and position oper- 
ators respectively) whereas that of the time operator T by jo — )■ —p, x ^ x,i ^ —i (for 



more details on VT symmetry see 16| and references there). The reason for choosing this 
particular potential (10. ip lies in the fact many physically interesting results jl|, S have been 
obtained while considering this potential in optical lattice. In particular, in ref [2|, by using 
spectral techniques and detailed numerical calculations the VT threshold Vq^ = 0.5 has 
been identified for this potential. It is found that below this threshold all the eigenvalues 
for every band and every Bloch wave number are real and also all the forbidden gaps are 
open whereas at the threshold Vq^ = 0.5, the spectrum is whole real line and there are no 



band gaps 17|. On the other hand, when Vq > 0.5 the first two bands (starting from the 
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lowest band) start to merge together and in doing so they form oval-hke structure with a 
related complex spectrum. The real as well as the imaginary parts of such a double valued 
band when Vq = 0.75 are depicted in Figure 2. We shall argue that one can anticipate by 
analytical argument that when Vq < .5 the VT symmetry of the potential is unbroken. We 
shall also see that there is another critical point (Vq) after which no part of the eigenvalues 
and the band structure remains real. 

Let us consider the following linear eigenvalue problem frequently used to study the linear 
properties of a periodic potential in optical lattice 

H^{x) = ^-j^ + Vix)tlj{x) = PiPix) (0.2) 

where (3 represents the propagation constant (energy eigenvalue) in the periodic structure 
and V{x) is given in (10. 1|) can be written as: 



V{x) 



2 + 2^1 - AVi cos {2x - i tanh-^2Vo) , Vq < .5 

2 + 2e2^^ Vo = .5 (0.3) 



2 + 2i\/AV^ -1 sin{2x-i tanh'^— ) , Vq > .5 



Recently in ref |l8l . Il9| , it has been shown that the concept of VT symmetry has its roots 
in the theory of pseudo-Hermitian operators . The reality of the spectrum is ensured if the 
Hamiltonian H is Hermitian with respect to a positive definite inner product {{il^i\ip2))r] '■= 
(?/'i|r7|'02), Vl'i/'i), |'02) € [a Hilbert space in which H is acting], where the Hermitian linear 
automorphism : "H — )■ "H is positive definite. Consequently, rj pseudo-Hermiticity of the 
Hamiltonian H is given by 

= rjHr]-^ (0.4) 



It has also been shown that 20| for any VT symmetric Hamiltonian H having real 
eigenvalues there exist an equivalent Hermitian Hamiltonian h by means of a similarity 
transformation h = pHp~^, where p = ^Jr\ is unitary [in the sense that it is Hermitian 
and for any pair of state vectors |0) : (('?/'|0))r? = ("01^0) = = (pV^Ip^) holds] 

which maps 'Kphys — ^ "H \Kphys being the Hilbert space equipped with the inner product 
defined above]. The Hermitian Hamiltonian h is equivalent to H and that they have same 
eigenvalues. 
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Now it is to be observed from fl0.2p and (10.31) that 
(i) for Vo < 0.5 the Hamiltonian defined in (10.21) satisfy the relations^ i/^ = 'qHrj'^ and 
h = pHp~^, where h = + 2 + 2a/1 — 4V^^ cos 2x, p = e~^^ [is a unitary operator in the 
sense stated earher], r] = p^ = e~^^ [is a positive definite operator^], 6 = tanh~^2Vo G M and 

Hence the Hamiltonian H is r^-pseudo-Hermitian with respect to a positive definite oper- 
ator rj = e~^^ and the Hamiltonian is linked with that of a Hermitian Hamiltonian h via an 
imaginary shift of coordinate which indicates that all the eigenvalues of H are real and the 
eigenvalue problem (10. 2p reduces to solving the Mathieu equation as shown below. 

(ii) For Vo > .5 we have rjV^x)!]'^ = V{x + iO) ^ V*{x), where rj = e~^^,9 = tanh^^^ 
which indicates that in this case the Hamiltonian H is not rj pseudo-Hermitian so the reality 
of the eigenvalues is not ensured. 

(iii) 



For Vo = .5, the potential V{x) =2 + 26^*^' had been studied in [ui ll7|, |22|. In 



ref. 17|, it was shown that the spectrum of this potential that is defined on the real line {x 



is real) is purely continuous and fills the semi axis [0,oo). There are no band gaps at all. 



This is consistent with the results of ref.[22||. In the latter paper it was also shown that the 
spectrum of the Schrodinger Hamiltonian for the above mentioned potential can be discrete 
if one chooses the complex contour appropriately. 

Hence it is clear that Vq = Vq^ = .5 is the transition point after which the eigenvalues of 
equation (10. 2 p becomes complex. This VT threshold Vq^ has been determined earlier in ref. 
Ill, l2| using numerical techniques. 

After making the change of variable z = x — ^tanh^^ ^ in case of Vq > .5, the equation 



(10. 2 p transforms into 

Under simple coordinate translation z ^ ^ — z the equation (10.51) reduces to Mathieu 
equation ^ + {a—2iqcos2z)ilj = 23], with characteristic value a = 2-/3 and q = ^/ AVq — 1. 



+ I2- 13 + 2i^AV^ - 1 sin2z] ^{z] /3, Vq) = 0. (0.5) 



A variety of numerical methods have been employed to solve Mathieu equation with an 
imaginary characteristic parameter as is the case in f|0.5D . Tabular values from an old paper 



of MulhoUand and Goldstein 2^ and Bouwkamp [25| refiect that the characteristic values 



q would be both real and conjugate complex: the transition occurring in the neighborhood 



^ In ref. 21|, it has been shown that the Hermitian hnear automorphism r] = e affects an imaginary shift 

of the co-ordinate: e~^'''xe^P = x + i9. 
^ ?7 is square of a Hermitian operator, so it is positive. 
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of g = go = 1.46876852 for solution having a period vr. The same critical (exceptional) point 
was obtained in ref.(26|, while considering the large limit of the quasi exactly solvable 
VT symmetric potential V{x) = —{i^sin2x + N)"^. So from equation flO.Sp it follows that 
there is a second critical point in the neighborhood of Vq = Vq ^ .888437 after which all the 
eigenvalues of the equation (10.51) becomes complex and occur in pairs as complex conjugates. 
Consequently no part of the band structure remain purely real. In the Figure 3 we have 



27 



plotted the band structure for the value of Vq = .8885 using Floquet analysis 

Now to determine the eigenvalues of the equation (10. 5p and correspondingly the band 
structure of the potential given in (10. ip . we proceed as follows |28|. As the coefficient of the 
eq. (l0.5p are vr periodic, so according to Floquet-Bloch theorem [27| the eigenfunctions are 
of the form t/j^z; /3, Vq) = (pui^] f3, Vo)e*'^^, where ^^(z + tt) = (p^iz) and u stands for the real 
Bloch momentum. Since <pu{z) is periodic so the entire solution of equation (10. 5p can be 
expanded in the Fourier series 



fcez 



which, inserted into eq. (10.51) . gives the recurrence relation 



C2k + C2fc[c2(fc-1) — C2(fc+1)] = 0, C; 



2k 



\/keZ 



(0.6) 



(0.7) 



[(z/ + 2A;)2 + - 2] ' 

for the coefficients C2fc(/3, z^). For any finite(truncated) upper limit n E Z, eq. (10. 7p can be 
written ClS Sb matrix equation 

Mn{iy;(3,Vo)Ci = (0.8) 
C2n) and A^„(z/; (3, Vq) is a 2r;, + 1 x 2n + 1 matrix given by 



where C„ 



[C-2nj 



■•) C_2, Co, C2, 

/ 1 

C2n-2 



Mn{iy;(3,Vo) 



-C: 



2n 



C2 1 -C2 

Co 1 
C-2 





-Co 

1 






-C-2 



v 



C- 



2n 



-C-2n+2 
1 



(0.9) 



For the finite dimensional case, finding the non-trivial solutions of the homogeneous system 
of eqn. (l0.8p is equivalent to solve 

det{Mn)=0. (0.10) 

The solution of the above eqn. (lO.lOp gives the relation between the energy eigenvalues /3 and 
bloch momentum i^. Figure 2 and Figure 3 have been drawn taking n = 3. 
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In conclusion, 

(i) we have shown that for a positive definite operator rj the potential (10 .ip is rj pseudo- 
Hermitian for Vq < .5 which ensures that all the band edges are real. It has been also that 
for Vo < .5 the eigenvalue problem (I0.2p reduces to solving the Mathieu equation. With 
the aid of rj pseudo-Hermiticity, the VT threshold for the potential under consideration, is 
obtained as V^*'^ = .5. This was previously obtained in ref. l2| using numerical technique. 

(ii) The existence of a second critical point ~ .888437 is established. Beyond this 
critical value no part of the band structure remains real. Numerical calculations also support 
this as is shown in Figure 3. At this point, it is important to note that in ref. 0], the authors 
have shown that it is possible to find stationary self trapped modes (with real eigenvalues) 
of the nonlinear equation ^ + V{x)il) + \'>p\^ip = (iip above the VT threshold Vq^ = .5. 
This is because of the fact that the band structure still remains real for some values of 
bloch momentum u even above the VT threshold. Stability analysis however revealed that 
the above class of lattice solitons are unstable. The possible physical implication of the 
second critical point Vq may be that there would not exist any soliton (neither stable nor 
unstable) solutions of the above nonlinear equation beyond Vq. It will be interesting to study 
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this physical aspect and other related physical properties of the optical lattice beyond this 
second critical point Vq. 
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